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Abstract. We extend the critical point self-consistency method used to solve field theories at 
their d-dimensional fixed point in the large N expansion to include superfields. As an application 
we compute the /3-function of the Wess-Zumino model with an O(N) symmetry to 0(1/N 2 ). 
This result is then used to study the effect the higher order corrections have on the radius of 
convergence of the 4-dimensional /^-function at this order in 1/N. The critical exponent relating 
to the wave function renormalization of the basic field is also computed to 0(1/N 2 ) and is shown 
to be the same as that for the corresponding field in the supersymmetric 0(N) a model in d- 
dimensions. We discuss how the non-renormalization theorem prevents the full critical point 
equivalence between both models. 
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1 Introduction. 



The Wess-Zumino model is the simplest four dimensional interacting quantum field theory which 
possesses supersymmetry, [jlj]. This bose-fermi symmetry is believed to play a role in the uni- 
fication of the forces of nature beyond the energy scales accessible by the present generation 
of accelerators. Nevertheless theories with supersymmetry are also of interest for testing out 
and developing new ideas. One reason for this is that one can calculate much more easily in 
the quantum super symmetric theory compared with say its bosonic sector, due to cancellations 
between graphs involving bosons and fermions. For example, there is no two loop term in the 
MS f3- function of the two dimensional super symmetric a model on an arbitary manifold, Q. In 
the corresponding bosonic a model the two loop contribution to the /^-function is non-zero, ||. 
Moreover, in certain situations, like the Wess-Zumino model, the renormalization is constrained 
in such a way that the calculation of the coupling constant renormalization is reduced to ob- 
taining the wave function renormalization. In other words ensuring unbroken supersymmetry in 
the quantum theory restricts the form of possible counterterms in the superpotential and hence 
determines a renormalization constant, |]l|, |4|. This non-renormalization theorem allows one to 
compute higher order corrections rather easily. For instance, the four loop /3-function for the 
pure model with one superfield is known, P, |, [7|], as well as that for its extension to include su- 
perfields with general internal symmetries, p|. Although the Wess-Zumino model was originally 
formulated in terms of ordinary spin-0 and | fields, there is a further calculational tool which 
substantially reduces the amount of work needed for these higher order calculations. This is at- 
tained by grouping the component fields in the same multiplet into superfields where the Lorentz 
space-time is extended to include Grassmann coordinates. Consequently in this superspace the 
Feynman supergraphs group classes of component field diagrams together thereby reducing the 
number that need to be computed. Moreover performing the integrations over the Grassmann 
variables, known as the D-algebra, the resulting Feynman integrals are invariably simple scalar 
integrals. This property can be better appreciated in superfield calculations in supersymmetric 
gauge theories where the number of component fields is larger than in the Wess-Zumino model 
and high loop integrals involving the component gauge field become extremely tedious compared 
with a scalar integral. 

Recently some insight into the structure of the Wess-Zumino model /^-function beyond the 
present four loop level has been gained in || through using the conventional large N bubble 
sum method. The leading order corrections in 1/N were computed for the case where the Wess- 
Zumino model has a multiplet of O(N) fields coupled to a separate scalar superfield. One of 
the motivations of that study was to ascertain the large coupling constant behaviour of the 
/^-function in a four dimensional model and hence to understand some of the problems in its 
coupling constant resummation, as well as being the foundation for repeating the exercise for a 
supersymmetric gauge theory. As leading order calculations do not reveal a substantial amount 
of structure it would be interesting to push the 1/N calculations to a higher order in both the 
Wess-Zumino and supersymmetric gauge theories. However, a formalism does exist to achieve 
this and is based on the critical renormalization group and (i-dimensional conformal field theory 



methods which were pioneered in the series of papers of ref. Jl0| , 11 1 for the O(N) a model in 
two dimensions. Before attempting to apply such methods to these models there are several 
problems which need to be solved. Although the methods of [JO], 11] have already been applied 
to the supersymmetric O(N) a model in |fl2[ [b| to obtain the /3-function at 0(1/N 2 ) and the 
supersymmetric CP(N) a model at 0(1/ N), [jl4j, these calculations were performed using the 
component lagrangians. Indeed the computation of 77 component Feynman diagrams in 
deduce the (3- function at 0(1/N 2 ) only serves to indicate that the development of a superfield 
approach would yield a more compact approach. This is one of the aims of this paper and will 
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be provided for the Wess-Zumino model. As an application we will then compute its /3-function 
at a new order, 0(1/N 2 ), as a function of the space-time dimension, d. Briefly, the method 
involves determining the critical exponents of the theory at the phase transition defined by 
the non-trivial zero of the cf-dimensional /3-function. As there is no mass at such a point the 
propagators of all the fields, for instance, involve the momentum raised to some power which is 
known as a critical exponent. From a field theoretic point of view this is a fundamental object 
as it is renormalization group invariant and therefore physically measurable. It depends only 
on the space-time dimension and the parameters of any internal symmetries. For our purpose 
this will be N. Therefore it can be expanded in powers of 1/N when N is large and deduced 
order by order by examining the scaling behaviour of the particular Schwinger Dyson equations 
truncated at the appropriate order. Moreover, as the exponents will depend on d = 4 — 2e they 
can be related through the critical renormalization group to the corresponding renormalization 
constants one ordinarily computes in perturbation theory. Hence one can deduce new coefficients 
in the perturbative series in successive orders in 1/N by performing the e-expansion of the 
exponents beyond the order currently known. 

Finally, another motivation for studying the Wess-Zumino model at this <i-dimensional crit- 
ical point is to understand the universality class to which it belongs. By this we mean the 
following. The bosonic sector of the O(N) Wess-Zumino model is effectively </> 4 theory with 
an O(N) symmetry. The field which is not an element of the multiplet O(N) acts as auxiliary 
which if eliminated yields the usual quartic interaction. As is well known (see, for example, 
fllql ) in ^-dimensions O(N) 4 theory lies in the same universality class as the lower dimensional 
O(N) bosonic a model and therefore the critical exponents one computes in either model are 
equivalent. Indeed these are of interest because both models underlie the physics of the phase 
transition of the Heisenberg ferromagnet. Therefore we are motivated to determine if there is 
an analogous lower dimensional equivalent theory for the Wess-Zumino model. Intuitively one 
would expect it to have some connection with the supersymmetric 0(N) a model. To achieve 
this we need to compare the values for the Wess-Zumino exponents with those of Jl^, [DjR . 
Such an equivalence would therefore play a role in some supersymmetric extension of the Heisen- 
berg ferromagnet in three dimensions, though of course it is not clear if such an object exists in 
the real physical world. 

The paper is organised as follows. We review the relevant properties of the Wess-Zumino 
model in section 2 and introduce the superfield approach for the critical point large TV" method. 
This is used to reproduce the known results at leading order in 1/N in section 3 where the 
0(1/N 2 ) correction is also derived for the wave function renormalization. The higher order 
correction to the /^-function is determined in section 4 where the technical details of the compu- 
tation of the three and four loop diagrams which occur are also discussed. Finally we collect our 
results in section 5 and produce several new coefficients that will appear in the MS /3-function 
at five and six loops. 



2 Critical point formalism. 

We will be working with an O(N) symmetric Wess-Zumino model containing one chiral superfield 
a and N chiral superfields 3> l , 1 < i < N. Its action is given by 



S 



d d x 



d 4 e + 



9 Z 



d 2 6o<$> 2 



d 2 9a$ 2 



(2.1) 



The coupling constant g has been rescaled into the cr-field so that the 3-vertex is in the form 



which easily allows us to use the technique of uniqueness, [17, [Tl[], to compute the higher order 
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graphs which will arise. Moreover the chirality of the fields in ( j2.1|) will prove to be useful in 
substantially reducing the number of possible diagrams which we will need to consider then. 
Before carrying out a dimensional analysis of the action to ascertain the canonical dimensions 
of the fields and the couplings, which is the first step in the critical point formalism, we recall 
some of the practicalities of the calculation of the perturbative /3-function. If we define a new 
field a 1 by a' = a/g then one would recover the conventional form of the interaction term. The 
non-renormalization theorem of ||| implies that the vertex is finite and so the overall vertex 
renormalization constant involves no infinities. Therefore the renormalization group functions 
of the rescaled interaction and its constituent fields satisfy the following simple relation 

(3(g) = [27* (p) + 7/ (<?)]<? (2-2) 

which implies the wave function renormalization purely determines how the coupling constant 
runs. The result (^^) was used in [||, [f| fjj to deduce (5(g) for the model without any internal 
symmetries. As we are dealing with the model with an 0(N) symmetry we require the pertur- 
bative results as a function of N. To achieve this we have specialised the four loop result of 
[]|] to the O(N) model. Useful in checking the relevant symmetry factors in this case was the 
package Qgraf, |fL8jl . Therefore, the MS results are 

7*(<?) = 9 ~ \{N + 2)g 2 - ~ (iV 2 - ION - 4 - 24C(3)) g 3 

-^(3iV 3 + 16N 2 + 152 N + 40) - ^-(N 2 - AN - 36) (N + 4) 



and 



3C(4)(JV + 4) + 20C(5)(7V + 2) ) g 4 + 0(g b ) (2.3) 



lA9) = y<? - Ng 2 + ^(2N + 1 + 6C(3)) 5 3 

+ N (j^iSN 2 ~ 56iV - 16) - ^C(3)(A r2 + 6iV + 44) 

+ ~C(4)(iV + 4) - 10C(5)(JV + 2)) g A + 0(g 5 ) (2.4) 

Hence, ( ^ ) implies 

(3(g) = \{d-A)g + \{N + A)g 2 - 2(N + l)g 3 + Q(iV 2 + 117V + 4) + 3(N + 4)C(3)) g 4 

+ Q(^ 3 - 36A^ 2 - 84iV - 20) - 3(iV 2 + 16iV + 24)C(3) 

+ + 4) 2 C(4) - 10(N + 2)(N + 4)C(5) \ g 5 + 0(g 6 ) (2.5) 

We note that for various reasons we have suppressed the usual factors of n associated with each 
power of the coupling in these expressions since, for example, they can readily be restored by a 
simple rescaling. Another reason is related to the critical point approach where we deduce critical 
exponents which correspond through the renormalization group to these renormalization group 
functions. The critical exponents are more fundamental quantities than the renormalization 
group functions in that they are renormalization group invariant. Therefore they take the same 
values in any mass independent renormalization scheme. This criterion arises from the fact 
that at a fixed point of the /3-function there is scaling and conformal invariance which would 
be broken by any non-zero mass. Therefore since, for instance, MS and MS are two examples 
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of mass independent renormalization schemes the critical exponents deduced in either scheme 
will be the same. The only difference is that the value of the critical couplings in each scheme 
will not be the same but they will be simply related by a constant rescaling of either of the 
couplings. However, throughout we will always have in mind that our results will relate to the 
MS scheme. 

With these renormalization group functions we can now develop the critical point formalism. 
We take the d-dimensional /3-function of (|2.5| ) and compute the location of the non-trivial d- 
dimensional fixed point as 

9c = | + (-& + !*> -8^-^ + 0^))^ 
32e - 176e 2 - 8(6C(3) - 37)e 3 



+ - (60C(5) - 9C(4) + 18C(3) - 4)e 4 + 0(e & ) j ^ + O j (2.6) 

in powers of 1/N where d = 4 — 2e. Although supersymmetry is not easily defined in d- 
dimensions we note that this is a valid step since in perturbative MS calculations with a di- 
mensional regulator ( p.5[ ) would be produced as the penultimate step in determining the 4- 
dimensional /3-function. Moreover it was shown in p2] , |I3| ] that the ci-dimensional version of 
the supersymmetric a model preserved certain features of the symmetry at its corresponding 
fixed point and this allows us to make similar assumptions about the Wess-Zumino model in 
(^-dimensions. For instance, the component fields of each boson and fermion partner in a su- 
permultiplet retained the same anomalous dimension in arbitrary dimensions. With (|2.6j) we 



can now dimensionally analyse the d-dimensional action (2.1) at g c choosing to work with the 
(T-field. Unlike the component calculation of [^] we define the canonical dimensions of the 
fields in relation to the full superspace and take the dimension of the Grassmann coordinates as 
— 1/2. As the renormalization group functions are non-zero we must also allow for the possibility 
of anomalous contributions which we will define with respect to the usual conventions of [ [Tcfl . 



Defining the full dimension of the <& superfield as a then the quadratic term of (2.1) implies we 
set 

a = /i - 1 + | (2.7) 

where r\ is the critical exponent corresponding to the wave function renormalization at criticality. 
In other words r] = y$(g c )/2. From the interaction term, if we denote by the full dimension 
of the a superfield and by x the vertex dimension, then 

P = 2^ - 1 - 2a - x (2-8) 

However, as the non-renormalization theorem implies the vertex is finite to all orders in pertur- 
bation theory then the (composite) operator a& 2 can have no anomalous piece. Therefore we 
set x = 0- This observation will have a simplifying effect in, for example, the computation of 
the 1/N 2 corrections to the Schwinger Dyson equations we will solve. Hence 

P = 1 - V (2-9) 

and the vertices of all the Feynman diagrams will in principle be one step away from uniqueness. 
From the remaining term of the action we deduce a scaling law relating the dimension of the a 2 
composite operator to the /3-function exponent. As P(g c ) = this is defined to be uj = — P'(g c )- 
Therefore, we have 

u = 2/i - 2 - 2(3 + r) a 2 (2.10) 



5 



where rj a 2 is the anomalous dimension of the bare composite operator. Although one can use 
( 2. ICQ to determine the /?- function by calculating the anomalous dimension of a 2 we have chosen 
to use the method of ]10| , 11] of solving the Dyson equations. It turns out that for scalar theories 
at least, both ways are related in that the Feynman diagrams where one inserts a 2 in a a 2- 
point function and examines the resulting diagrams, it is easy to see that they have a complete 
analogue with the diagrams which arise in the Schwinger Dyson approach. However, both 
methods require a similar amount of effort to determine the solution. 

Near g c , we can define the asymptotic behaviour of the dressed propagators of and a in 
Euclidean space. As p — > oo, 



(*(-p,0)*(p,0O) 



A5\e-e') 

(p 2 y- a 

B5 A {9 - 9') 

( P 2 y-p 



(2.11) 



where we have set d = 2fi and A and B are unknown momentum independent amplitudes. 
Corrections to these scaling forms can also be defined. For example, we take 



{a(-p,9)a(p,9')) 



A5 4 (9-9') 
B5\9 - 9') 



1 + 



1 + 



A 1 



{p 2 Y 

B' 



{p 2 Y 



(2.12) 



In principle one can include other correction terms involving other exponents with different 
canonical dimensions. However, renormalizability and hyperscaling laws ensure that knowledge 
of two independent exponents such as r\ and v is sufficient to deduce the others. Here uj is 
the exponent which relates to corrections to scaling. The two pairs of amplitudes (A, B) and 
(A 1 ,B') are independent of one another. Following the methods of 1C, 11 1, we will use ( |2.12| ) to 
solve the Dyson equations. The inclusion of a non-zero anomalous dimension in the exponents 
of the propagators means one considers only those diagrams with no dressing on the propagators 
which of course will further reduce the number of graphs which would need to be computed at 
high orders. 

Ordinarily in the critical point method of |T(], [TT| t he next stage is to determine the asymp- 
totic scaling forms of the inverses to ( 2.11 ) and ( |2.12D by inverting them in momentum space. 
This is because they will appear in the Dyson equation as illustrated in fig. 1. However, the 
situation is different here and can be best introduced by first recalling the nature of perturbative 
calculations in superspace. In the calculation of the corrections to, say, the quadratic terms in 
the effective action one needs to retain the external momentum space superfields explicitly when 
representing the graph from the super-Feynman rules. The reason for this is that the interactions 
and propagators involve supercovariant derivatives which act on the Grassmann coordinates and 
hence the external fields since they too depend on 9. In ordinary non-supersymmetric calcu- 
lations it is not necessary to include the external fields since it is the case that for (local) 
quantum field theories derivative couplings are defined in coordinate space. Therefore when 
performing computations in momentum space their Fourier transform corresponds to vectors of 
the conjugate variable which therefore do not act differentially on the external legs. After the 
manipulation of the supercovariant derivatives, known as performing the -D-algebra, one is left 
with a simple momentum space Feynman integral which can in principle be computed. The final 
result is an integral over the full momentum superspace which includes the fields corresponding 
to the external legs. For a 2-point function this would be of the form, p9l |(|, 



d d p 



d 4 9f( P M-p,9Mp,9) 



(2.13) 
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where <3? and $ are the external fields whose momentum is p. Therefore this would be regarded 
as the correction to the corresponding term in the effective action. From the point of view 
of perturbative renormalization the extraction of the pole part in f(p) with respect to the 
regularization would then be absorbed by the appropriate counterterm. In light of this the 
critical point approach to the Dyson equations is adapted to be similar. One manipulates the 
graphs in momentum superspace using ordinary D-algebra but with the lines of the graphs of fig. 
1 replaced by the critical propagators, Q2.ll ) and ( |2.12 ), and the external fields assumed to be 
included. One again obtains an expression of the form (|2.1^ ) but the corresponding f(p) will now 
depend on the exponents of the lines, a and /3, of the original graph. The momentum dependence 
will be in scaling form with the exponent given by the dimension of the original integral. The 
remaining part of ( 2.13| ) would be regarded as the propagator. In the approach of [10, 11] in 
non-supersymmetric calculations this would not appear here since the fields of the external legs 
are not included which is the reason why one needs to invert the corresponding (2.11). So in 
the superspace approach this is not necessary since the propagator appears naturally with the 
scaling part of the graph. Although we have included the inverses in fig. 1, these are also meant 
to account for the situation when one has the corrections to scaling of ( p. 12 ). Therefore, we 
define 



3 _1 (p) 



1 

A 



1 



A' 



(p 2 



a 



-i 



to 



1 

B 



1 



B' 



(P 2 



(2.14) 



which is similar to the situation in |l0| , O]. In practical calculations it will be a simple exercise 
to verify that the above argument leads to dimensionally consistent equations. 



3 Solving the Dyson equations. 

We are now in a position to determine expressions for the critical exponents and begin with 
the leading order analysis. With the propagators from ( 2.12] ), the Dyson equations, which are 
given in fig. 1 and have been truncated at leading order in 1/N, may be written in the following 
integral form in the scaling region as 








cr- l {p) + ^NA 2 H(a,p) J d 2 ^k 



1 + 



A' 



+ 



B' 



{k 2 y- a {{k - p) 2 y- 



i + 



A' 

Jk^y 



{{k-p) 2 y 

A' 



{{k- P ) 2 y 



(3.1) 



where the functions Ti of the fields and external momentum represent the 6 integral of ( |2.13| ) 
and we have omitted terms which are quadratic in A' and B' . We note that we use a form of the 
Dyson equations where the propagators are dressed but not the vertices. The factor of § comes 
from the symmetry of the first graph in fig. 1. Substituting for a, f3 and uj and decoupling (|3.l| ) 
into terms with the same momentum dimension, we obtain 








1 + 
1 + 



zv (fj, 

\Nzv (2fi 



1 + r], 1 - §t?) 
2 + 77, 1 - l v , 1 - ±77) 



(3.2) 



1 + J - \r],n 



[u (fl 

+ V (2/i — 3 + id' + 7], 2 — LU' 
zNv (yU — 1 + U) — |r7, 1 — 577, fl + 7] — lo') A' 



1 + 77, 2 - i?y - u') z - 1] A' 
±77, 1 - ±7?) zB' 



B' 



(3.3) 
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where we have set z = A 2 B and defined the anomalous piece of uj by uj = /i — 2 + uj' . The function 
1/ is defined for arbitrary a, j3 and 7 as z/(a,/3, 7) = a(a)a(f3)a('y) with a(a) = T(fi — a)/T(a). 
Eliminating z from (|3.2[) we obtain the relation that will give us 77 at leading order in 1/N 



4T(1+t//2) 
NT(ja-ti/2) 



v(2- (i-ri,n 



1 + ±77, /j 



1 + 7?) 



Setting 77 = and expanding (|3.4|) in powers of 1/iV we find 771 is 



m 



4T(2fi - 2) 



r( M )r2(Ai-i)r(2- M ) 



(3.4) 



(3.5) 



Remarkably the computation of 7/2 requires only the expansion of ( |3.4[ ) to the next order in 
1/N. Ordinarily one would have two and three loop corrections to each of the Dyson equations 
of fig. 1 but in the Wess-Zumino model the chirality condition on the vertex excludes these. 
Therefore the first non-zero corrections will arise at 0(l/iV 3 ). In this respect the Hartree-Fock 
approximation at 0(1/ N 2 ) determines the full value of 772. Therefore we simply obtain 



m 



+ ^(2/i-2) - ^(M-l) - ^(1) + 



1 



2(/i - 1) 



(3.6) 



where ip(x) is the logarithmic derivative of the Euler T-function. Moreover the absence of these 



higher order contributions means that unlike the bosonic a model calculation of [10, 11 1 the 
Dyson equation is finite to this order and does not need to be either regularized or renormalized. 
Regularization is usually necessary due to vertex subgraph divergences in the higher order graphs 
as the vertices of the Dyson equation are undressed. It is introduced by shifting the analogous 
o"-field exponent by an infinitesimal amount, A, which is equivalent to setting \ ~^ X + A. 
As there is no vertex renormalization in the Wess-Zumino model and since the critical large N 
method is effectively a perturbative expansion in the vertex anomalous dimension then we are 
led to reason that at the fixed point the Dyson equations do not require any renormalization at 
any order in 1/N. For subsequent calculations we will require the value of z to two orders. Also 
setting z = J2^=i Z i/N l knowledge of 771 and 772 allows us to deduce 



Zl 
Z2 



- 5 r (M)m 

\t<jjl) m 



li) + V(2M-2) - V^-l) - ^(1)]t7? 



(3.7) 



determines corrections to the exponent uj. If one formally 
= (3.8) 



The second pair of equations, 
writes ( |3.3j ) in matrix form 



then uj is determined by requiring that the determinant of this matrix vanishes, [11]. This gives 




1 = z 2 Nu (2/i - 3 + uj' + 77, 1 - ±77, 2 



UJ 



§77) v (fi - 1 + J - ±77, 1 - §77, fi + 77 



UJ 



+ ZV (/i — 1 + UJ 



1 + 77,2 



±77 



UJ 



(3.9) 



However, in deducing ( |3.S| ) one needs to bear in mind that the extra graphs at 0(1/N 2 ), which 
would usually play a role in the evaluation of 772 but which are absent due to chirality, would 
also in principle need to be included here. Though they would only enter in one of the Dyson 
equations because the leading iV-dependence of each of the entries of the matrix of ( |3.8|) , when 
considered in coordinate space, is not the same. Moreover this feature has already been observed 
and discussed in similar calculations in other models, ]^0|. Although the relevant graphs are 
absent in determining the first non-trivial term of uj certain topologies which cannot be ruled out 



S 



by chirality occur at higher order and therefore will need to be included in the relevant Dyson 
equation of ( |3.8| ). Therefore setting u>' = YliLi u i/N l and retaining only the leading terms in 
(|3,9|), we obtain 

_ (2\x - 3)771 _ (/i - 1)771 

(wi-r/i) 2( M -2)iV 1 ' J 

As the second term is suppressed by a factor of 1/N compared with the first we can ignore it 
and rearrange ( 3 . 1 0| ) to find 

ui = - 2(m- 2)771 (3.11) 

We need to compare these results for rji and u% with the results already obtained via the 
conventional bubble sum methods, ||, and demonstrate their equivalence. To do so we must 
convert the critical exponents into the original renormalization group functions. From the defi- 
nition of the critical coupling, (|2.6|) , we can replace the /j-dependence in (|3.5|) at leading order 
by setting fi = 2 — Ng c /2 which determines 7<j>,i(fl , c ) as an explicit function of g c since 7$,i(<7 c ) 
= rji/2. Then it is easy to deduce for non-critical couplings that 

. 2r(2 - y) ( 1 \ 

= ivr(2- y /2)r 2 (i- y /2)r( y /2) + °{w) (3 ' 12) 

where we have introduced the scaled coupling y = Ng for convenience. For subsequent calcula- 
tions it is also useful to define a new function G(y/2) by 

G (y_\ = r ( 2 ~^) f3131 

\2] r(2-y/2)r 2 (i-y/2)r(i + 2/ /2) 1 ■ ] 



Taking into account the symmetry factor of 1/2 present in each $ loop, (3.12) is precisely the 
result obtained in ||. We can perform a similar exercise for the /3-function. Noting that uj\ = 
2y 2 c G(y c /2) and defining (i{g) at 0(1/N) as 

P(g) = (//-% + \(N + A)g 2 + 5 2 /i(^) (3.14) 

where fi(gN/2) is an arbitrary function we see that at the critical point we must have 

= 2(/i-2) + (iV + 4) 5c + 2g c f 1 (^j (3.15) 

Using ( |3.15| ) to eliminate the (N + A)g c term in (3'{g c ) we find that 

m = - ^ (3.16) 



Therefore we find in four dimensions 



1 r9 iy / 2 / 1 \ 

(3(g) = -(N + A)g 2 - Ag 2 J dxG(x) + O f^J (3.17) 

which was obtained originally in Q. To compare with perturbation theory one computes the 
integral by first expanding the integrand in a Taylor series as a function of x before integrating 
each term individually. It is elementary and reassuring to check that the coefficients of (|2.3| ) 
and ( |2.5| ) are reproduced at four loops respectively from (3.12) and ( 3.17] ). It is worthwhile 



comparing some of the features of both methods that lead to the same results. In |9[ an infinite 
number of diagrams had to be summed which was only made possible by a seemingly miraculous 
cancellation of terms. By contrast one needs only to calculate two one-loop diagrams in the 
critical approach. The elimination of the amplitude factor z between both Dyson equations in 
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effect reproduces the conventional bubble sum giving the same results in the end. However, 
the advantage of the critical method lies in performing 0(1/N 2 ) calculations. Indeed it is not 
always clear if the corresponding bubble sums at this order would prove to be easily calculable 
but would certainly rely on further remarkable cancellations. We can continue the exercise 
of converting the exponents into the renormalization group functions by considering 772 • By 
retaining the higher orders in the critical coupling we find that to 0(1/N 2 ) 



f \ y nfy\ , 2 y 



nf) + M5)M§ 

+ yG 2 ( y -)(*(y) 



1 



(3.18) 



.27 V y " y (y-i) (y-2) 

where, for later convenience, we have introduced the function (5\(gN/2) which is defined by 

(3(g) = Gu-% + + g 2 (3i( 9 4-) +o(±) (3.19) 



2 * ' V 2 J \N 2 

and represents all the 0(1/N) corrections. Its explicit form can be deduced from ( |3.17| ). The 
function ^f(y) is given by 

*(y) = ^(1 - y) + V (l + I) - V (l " I) - V^(l) (3.20) 

and its expansion in powers of y will only involve the Riemann zeta series. The higher order 
correction to (|3.19|) will be determined later. If we expand (3.18) to 0(g 4 ) we recover all the 



coefficients of ( |2.3[) at 0(1/ N 2 ) which leads us to believe the result for 772 is correct. 

We can now address one of the issues raised in sect. 1 concerning the equivalence of the 
Wess-Zumino model with a lower dimensional theory. The criterion for determining which model 
this ought to be is dictated by the nature of symmetries and the interaction. At criticality the 
quadratic terms in an action primarily determine the canonical dimensions of the fields and 
couplings. Therefore in the case of the O(N) Wess-Zumino model an obvious candidate is the 
O(N) supersymmetric a model which is defined in two dimensions and has the superspace action 



s = - ~ / <r'.r <fo 



-D$D§ - a ( $ 2 - - 
2 V A 



(3.21) 



Here the coupling constant A differs from that of fl2.1|) in that it is dimensionless in two dimen- 
sions and its critical value is deduced from its /3-function which will clearly be different from 
(|2.5| ). As this model has also been studied to 0(1/N 2 ) we are in a position to compare the 
critical exponents of the fields in both models. (Comparing the exponents v = — (3'(X C ) of the 
supersymmetric a model with u = — (3'(g c ) of the Wess-Zumino model is not sensible as their 
respective canonical dimensions are (/J — 1) and (fx — 2).) From the results of []l(], ^] it transpires 
that the $-superfield dimension has the same (f-dimensional value to 0(1/N 2 ). Whilst this was 
expected at leading order because the diagrams are the same, the agreement at next order is 
remarkable when we compare the nature of the calculation of 772 in both models. In the Wess- 
Zumino model 772 is determined purely by iterating the equation which produces r/i to the next 
order because there are no new contributing graphs at 0(1/N 2 ). By contrast in the a model 
to obtain 772 one has to compute the contribution from two and three loop topologies which are 
divergent and need to be regularized in the manner discussed earlier. In the component field 



calculation this equated to evaluating 24 graphs, [12|. Therefore it would appear that both (|2J 



and ( 3 . 2 1| ) are equivalent at least at 0(\/N 2 ). In other words 7$(<? c ) = 7$(A C ) to this order. 



However, the cr-field of both models does not have the same dimension as a direct consequence of 
the non-renormalization theorem. As noted earlier x = in the Wess-Zumino model but in the a 
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model x is non-zero. For example, at leading order xi = — This property can be illustrated 
further if we consider calculating the exponent v = fi — 1 + 0(1/N) in the Wess-Zumino model 
and compare it with the value in the a model, [13]. To do this we can replace the exponent 
uj in the correction term in each of the scaling forms fl2,12|) with v. Consequently the leading 



order value v\ can be deduced from the equation analogous to ( |3.9| ) where now both terms on 
the right side will contribute giving 

m 1 (3.22) 



(r?i-2i/ x ) 



which implies v\ = 0. This is the same as the value in the a model, and would again 

have encouraged us to believe in a possible equivalence. However, unlike the situation with u 
the coordinate space version of the matrix consistency equation which determines v does not 
have the reordering problem and we can therefore also compute vi- This is essentially trivial 
due to the absence of the higher order two and three loop topologies due to chirality which 
again means there are no new higher order graphs which need to be included. Therefore v<i 
is deduced by expanding the consistency equation which gave v\ to the next order in 1/N, in 
the same way that 772 was deduced. We found that vi = 0. In the a model v has a non- 
zero contribution at 0(1/N 2 ), |13[| . It is not clear whether v will vanish to all orders in 1/N 
in the Wess-Zumino model since there will be new graphs at next order. Their contributions 
would have to cancel to have 1/3 = 0. In light of these remarks it would certainly be very 
surprising if 773 calculated in both models was the same since the effect of differing values of x 
would almost certainly be important at 0(l/iV 3 ). In all these observations concerning a possible 
equivalence the underlying reason for its failure after 0(1/ N 2 ) is due to the chirality condition 
which implies x = 0- However, another way of understanding it is by considering supersymmetry 
in relation to the space time dimension. Clearly the Wess-Zumino model is invariant under AT 
= 1 supersymmetry in four dimensions. However, theories in two dimensions which are built 
out of chiral super fields are invariant under J\f = 2 supersymmetry. Therefore if one wished 
to construct the lower dimensional equivalent theory to ( |2.1| ) it would need to have more than 
one supersymmetry. The 0(N) supersymmetric a model unfortunately does not satisfy this 
criterion being invariant under M = 1 supersymmetry which is why it is remarkable that its 772 
and that of the Wess-Zumino model are the same. 



4 Calculation of 0J2. 



We are now in a position to extend the result for u% to the next order in 1/N. Unlike the 
calculation for 772 which simply involved the expansion of ( [3.4D to next order the reordering 
which occurs in ( |3.3|) means we need to include the contributions from several higher order 
graphs. These are illustrated in fig. 2. In the analogous calculation of the j3- function in (f) 
theory, pi| , 22|, there were 33 graphs which needed to be computed. In principle the same 
graphs would have to be considered here but the chirality condition excludes the majority of 
these to leave the topologies of fig. 2. Consequently if we denote their total value by n2, the 
solution of the Dyson equation, (|3.9[) , can be expanded to the next order in 1/N to obtain 



u>2 



40i-l)0i-2)¥( M ) 



(16/j 4 - 96/u 3 + 188/j 2 - 136/j + 25) 
2(2/x- 3)^-1) 



+ n 2 



(4.1) 



where ^(/i) = i(^(2fi — 3) + ifj(3 — fi) — tp(^ — 1) — V'(l)- The terms not involving U2 correspond 
to the expansion of the a-functions of ( p.9[ ). We note that at 0(1/ N 2 ) we can no longer neglect 
the second term of ( |3.10| ). Before we can substitute the lines with the critical propagators, 
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we need to compute the D-algebra of each graph. Elementary application of the rules of || 
means that we are left with a D-operator acting on several lines which in momentum space will 
correspond to an extra factor of p 2 where p is the momentum flowing in that line. We note that 
in perturbation theory this invariably means that that line vanishes from the Feynman integral 
since massless propagators have unit exponent. In the critical large N case this will not always 
be the case as the D-operator can act on a cr-line whose canonical exponent is not unity. Before 
detailing the evaluation of each graph we record how the D-operators are distributed around 
each of the graphs of fig. 2. For the non-planar graph each of the upper $ lines joining the 
external a lines has a D-operator, which means in the light of the above remarks that they now 
have zero exponent. The analogous lines of the four loop graph also have one of these operators 
acting on them to also reduce them to zero exponent in momentum space, as does the completely 
central a line. 

It turns out that for the non-planar graph this D-algebra immediately reduces it to the simple 
two loop topology given in fig. 3 where we use the coordinate space convention of integrating 



over vertices in the language of [11]. This graph can be reduced to an integral whose value 



is well known if one applies the transformation y of [11]. The result is an integral originally 



computed by Chetyrkin and Tkachov, 23], denoted by ChT(l, 1) in the notation of [11]. Hence 
the full contribution of the non-planar graph to H2 is 

- |0*-1)6(m)i^ (4.2) 

where Q(fi) = ip'([J> — 1) — We note that each term of the e-expansion of 0(/i), with fi = 

2 — e, is proportional to the Riemann ("-function beginning with C(3). 

The treatment of the remaining graphs is more involved since there are three <r-lines where 
the correction in the asymptotic scaling functions can be inserted and the distribution of the D's 
around the graph introduce an asymmetry. Otherwise two of the insertions would be equivalent 
under a symmetry transformation of the integral. Therefore we cannot assume a priori that 
their contributing values would be the same. Of the three possibilities it transpires that that 
graph where the cr-line insertion is on the bottom line of the graph in fig. 2 (with the D's on the 
upper <£-lines) has already been computed in the bosonic 0(N) a model calculation of its two 
dimensional /3-function, plj| . For completeness, we note that the elementary rules of integration 
by uniqueness, [ p]] , after transforming to the momentum representation leaves the three loop 
integral illustrated in the first graph of fig. 4. Therefore, we merely quote the contribution of 
this graph is 



Qi-l)(2/j-3) 
2(/x - 2) 



Vt (4-3) 



where $(//) = ip'{2^ - 3) - ^'(3 — (i) — - 1) + ^'OO- Lik e (@3) th e e-expansion of the 
function of \x within the brackets also only involves the Riemann £ series and no rationals. 
However, here the first term involves C(5)- 

It turns out in fact that the graph where the cr-insertion is on the top cr-line relative to the 



□ 's is equivalent in value to (4.3), at this order in 1/N. This can be established either by explicit 
evaluation or by a series of conformal manipulations on the three loop integral which results after 
the same elementary steps which produced the first graph of fig. 4. Instead the second graph of 
fig. 4 is obtained. To relate it to the former we first perform a conformal transformation with 
the origin on the right which produces an integral with the same topology but with different 



exponents on the lines. Replacing the (y — u) line by a chain integral, [11], and choosing the 



exponents of these two lines in such a way that the top right internal vertex is now unique, one 
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obtains a three loop graph with the same topology but now with exponents (/i — 1) on all the 
lines except for the (y — u) and (x — z) lines whose exponents are (3 — n). Finally, by attaching a 
propagator of exponent (2/j, — 2) to the x external point using the chain integration rule of [|ll]], 
the final integration over this unique vertex yields the first graph of fig. 4. Moreover the factors 
associated with each of these integration steps reduce to unity. Therefore both these insertions 
are equal in value. 

The remaining integral with the insertion on the central cr-line cannot immediately be reduced 
to one whose computation has been given previously. Therefore we will give its evaluation in 
detail. Transforming to the momentum representation yields the first graph of fig. 5. After 
carrying out a conformal left transformation on it the resulting integral has a unique triangle, 
so that when it is completed the subsequent integral has a unique vertex. Performing this and a 
simple chain integral, before undoing the original conformal left transformation, one is left with 
the second integral of fig. 5. The associated factor is a 4 (/x — l)/a(2/j — 4) relative to the first 
graph of fig. 5. Then applying the successive transformations — > and <— in the language of |p4j| , 
before applying the momentum representation transformation yields the final graph of fig. 5. 
Ordinarily in perturbation theory when such a topology is encountered where all the lines have 
unit exponent, the way to proceed is to integrate by parts. Although the same strategy applies 
here, it turns out that the set of graphs which result from this are individually divergent, though 
the overall result is clearly finite. To handle this one has to introduce a temporary (analytic) 
regulator 5 in the power of various lines. The final result will of course be independent of 
any choice, but it is judicious to choose the distribution of <5's in such a way that subsequent 
integrations can be carried out. To this end we replace the exponents of the lines with (2/i — 3) 
by (2/j — 3 — 5) and the central lines with exponent 1 by (1 + 5). This symmetric distribution 
means that two integrals result after integration by parts, each of which can be immediately 
reduced to a two loop integral where there is a common factor proportional to 1 /5. After several 
elementary transformations each of these has the 5 = form ChT(3 — fj,,fi — 1). However, the 
1/5 pole means that each 5^0 integral needs to be determined to the 0(5) term. Unfortunately 
this is not possible for each case. Instead it turns out that due to the distribution of <5's in each 
integral we can compute the difference in the values of the integrals to the order we require and 
it is merely given by the value of ChT(3 — fj,, fj, — 1 — 5) to 0(5). This is achieved by comparing 
the corresponding exponents in each graph. In other words if the lines have the same exponent 
for non-zero 5, then to 0(5) they will not influence the value of the difference. Therefore in these 
lines we can effectively set 5 = 0. Only for those lines where there is a mismatch in the exponent 
in corresponding lines must one not set 5 = in them. After matching in this way, one essentially 
reduces the calculation to the evaluation of [ChT(3 — fx, \i — 1 — 6) — ChT(3 — /i, /x — l)]/5 as 5 
—* 0. Thus piecing the steps in this summary together, the total contribution to H2 from this 
integral is, including the amplitude factors, 



(/i-l)(2/x-3) 
2(n - 2) 



i, , i 2/ \ 8(2/i-3) 
$(/j) + tf 2 (/z) 



2)771 (ii-2f 



rfi (4.4) 



Although the function within the brackets is different from Q4.3| ) its e-expansion is similar in 
that it begins with C(5)- 

Finally, it is an elementary exercise now to substitute the values for these integrals into (4.1 
to find 



W 2 = Tfl 



3(3/x-4)(//-l)e(//) _ (2/i-3)(/i-l)fr 2 (/i) 4(2/1 - 3) 2 
2(/i-2) (m-2) ( M - 2)2^ 

(4/i 3 - 24/i 2 + 50/j - 35)Qu - 1)9 (fi) _ (2ji - 3)(ji - !)$(//) 
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+ 



(16/x 7 - 192// + 956// 5 - 2536// + 3825// - 3212/x 2 + 1328// - 182) 
20i-l)(/z-2)3(2//-3) 



(4.5) 



We draw attention to the fact that as in other models where there is a reordering of the Dyson 
equation to produce the /3-function, [20], there is a term linear in rji which prevents the 0(1/N 2 ) 
correction being proportional to rj 2 . 



5 Results and discussion. 



Following the same procedure we used to determine 7$(<?) to 0(1/N 2 ), we can write down the 
/3-function in a similar way. If we define the 0(1/ N 2 ) /3-function by 



(3(g) = (ji-2)g + 



Ng 



then a similar derivation to that which produced ( |3,17| ) gives 



2/3? 



1 rgN 



where 



(y-2)(y-l)y \^ 2 (y) + $(y) 

y( 4y 2 _ 7y + 4) y2 



l*(y) 



y + (y - 1) 



2ft 2 (y) 



l 



(5.1) 



(5.2) 



V 



(y-2) (y-i) 

and the new functions $(y) and G(y) are given by 

<&(y) = V'(i - y) - $ (i + 1 



+ ( 2 /-2)( ? / 3 + y 2 -6) 



G 



(y-i) 
2fy 



6(y) 



(5.3) 



^ (l - 1) + 1^(1) 



e(y) 



1^(1) 



(5.4) 



After expanding the integrand in powers of y and completing the integral, it is easy to verify 
that (|5,3j ) reproduces the first four coefficients of ( [2.5| ) at the appropriate order in 1/N which 
gives us confidence that ( [4.5[ ) is correct. Having therefore verified that our result is consistent 
with all known information we can now provide some new coefficients in the perturbative series 
in four dimensions beyond ( |2.5| ). For instance, if we represent the unknown 0(1/N 2 ) higher 
order contributions to the /3-function by 



then we can deduce 



° 5 = 16 



1 1 00 

-(d-4)g + -(N + 4)g 2 + ^(a r N + b r )N r - 2 g r + 1 

r=2 



C(3) 



(5.5) 



„ | C(3) 3C(4) 

6 40 20 40 



14 



^!+(^-»«3>-> + >-> 

6 ° = I" (^-I5« 3 » + >») C(3) + S C(4) 

57 1 53 

- yC(5) - — C(6) - HC(7) (5.6) 

The five loop coefficients, which are new, will provide a useful crosscheck on future full pertur- 
bative calculations in MS in this model. 

Ordinarily once an exponent has been checked for consistency with perturbation theory one 
can gain additional information about the same exponent in the three dimensional theory since 
the expression is usually valid in d-dimensions. However, attempting this for Q4.5| ) one finds that 
U02 diverges as fi — > 3/2. This is due to non-cancelling singularities arising from the functions 
$(/i), ^ 2 (//) and the final term. If one compares this to the situation in 0(N) </> 4 theory this 
is somewhat unexpected since its u>2 is well behaved in three dimensions, |2l], For ( |4.5| ) 
one can trace the source of the terms which are divergent in three dimensions to the original 
Feynman integrals. It turns out that they arise solely from the four loop integral of fig. 2. In 
the calculation of [21] there is an additional singularity from the graph which is a non-planar 



version of the four loop graph of fig. 2. Its contribution in O(N) c/> 4 theory precisely cancels the 
singularity from the planar four loop graph. As we have already noted the non-planar topology 
does not arise in the Wess-Zumino model due to chirality and it is its absence here which leads 
to a singular 14%. In making this comparison it should be pointed out that whilst the explicit 
/i-dependent values of the analogous four loop integral of fig. 2 will in general be different in 
both models, each will have the same combination of divergent functions, $(//) and ^ 2 (/z), with 
the same residues as \i — ► 3/2. 

In producing expressions for the ^-function like (|T| and Q we have so far assumed that 
they are always evaluated by first expanding the integrand in a power series about the lower end 
of the integration range. Each term of the series is then integrated separately. However, we can 
also regard these expressions as being a summation of part of the full (unknown) perturbative 
series for /3(g). Although the explicit integration cannot be performed even for /3\(gN/2) we can 
at least determine some information on the range of validity of these expressions and how the 
subsequent corrections affect it. As was pointed out in || the obstruction to the resummation 
at 0(1/N) is given by the first singularity in the integrand of (3\(gN/2). This was found to be 
at g = 3/N due to the singularity in the numerator T-function of G(y/2), fl3.13|) . Although this 
is also divergent at y = 2 there is a compensating pole in the denominator of G(y/2) at the 
same point. In other words, 

Therefore to 0(1/N) /3(g) diverges logarithmically to + oo as g — > 3/N. Moreover, it is also 
positive in this interval. Repeating this analysis now when the explicit expression for (32(gN/2) 
is taken into account we observe that the first singularity is now at g = 1/N It arises purely 
from the ^ 2 (y) and <&(y) terms of ^(y) and in particular 

1 98 1 

which again leads to a logarithmic singularity in (3(g) but now at g = 1/N and tending to negative 
infinity. Therefore the effect of including the higher order corrections is that the obstruction 
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to the resummation moves closer to the origin. Moreover, this singularity is also related to the 
absence of a finite value for 0J2 in the three dimensional model. This is simply because the leading 
order relation between the coupling and d is fi = 2 — y/2. It is also interesting to compare these 
results with the analogous ones for the 0(N) </> 4 /^-function at 0{l/N 2 ). From the expressions 
for loi and L02 given in [21, 22] and the perturbative /^-function of pij, we find that 0i has a pole 



for y = 15, whereas fa is singular for y = 9. Again, the pole in fa stems from the (analogous) 
\^ 2 (y) + &(y)] term. The reason why the coupling range is so much bigger in the (ft 4 model 
is simply that the coefficients of the one-loop terms (which determine the arguments of the 
functions G, Vl/, <1> and O) are 1/2 and 1/6 in the Wess-Zumino and (ft /3-functions respectively. 



Moreover the signs of corresponding integrand poles are the same as (5.7) and (p.q). Therefore 
in both models the radius of convergence of the large N /3-function decreases when higher order 
corrections are included. Further, at 0(1/ N 2 ) a spurious fixed point emerges since both models 
are non-asymptotically free and hence their /^-functions are positive for very small couplings. 
In explicit perturbative calculations similar zeros can also occur when higher order terms are 
successively added. Since their appearance fluctuates with the number of terms included in the 
series we would expect that the zero at 0(l/N 2 ) would disappear if the 0(1/N 3 ) correction was 
included. In this case it also would be interesting to ascertain where the first integrand pole 
occurs. 

We conclude with several comments. Although as we have discussed there is not a full 
equivalence between the O(N) Wess-Zumino model and the O(N) supersymmetric a model we 
can make some observations on the numerology of the wave function exponent in both models 
and examine the effect supersymmetry has on 772- By this we mean the presence of rational 
and irrational coefficients that appear in the corresponding renormalization group function. For 
example, the /3-function of the Wess-Zumino model clearly contains rationals and the Riemann 
zeta series. Likewise if one expands 772 to higher order rational numbers will occur. Although 
this may not appear surprising, the e-expansion of 772 for the 0(N) supersymmetric a model 
near d = 2 — e dimensions, which is the same as the result for the Wess-Zumino model, yields 
only the series £( n )> n > 3, beyond two loops to all orders in perturbation theory. The (5- 
function at 0(1 /N 2 ) has the same property, |13[ |. By contrast the bosonic O(N) a model 
has rational coefficients at two and higher loops in its renormalization group functions, || ||], 
and the imposition of supersymmetry renders the cancellation of the rationals by fermionic 
partner graphs. Naively one might expect a similar property for four dimensions. However, in 
rewriting the tft(x) functions of ( |3.6| ) to ensure their e-expansion has only C( n ) terms now with 
respect to d = 4 — 2e, one introduces functions whose expansion will produce rationals. This 
arises from the iteration of the Hartree Fock approximation and, at least for the Wess-Zumino 
/5-function, we have shown at 0(1/ N 2 ) that the rational coefficients in the four dimensional 
result also arise from a similar iteration. The higher order graphs included here involve only 
C(n)'s in their expansion as we emphasised in their computation. Indeed in light of these 
remarks several questions arise as to whether the imposition of another supersymmetry on this 
model would reduce the renormalization group functions of that model to rational at one loop 
only and irrational thereafter. Also, it is natural to question whether all the higher order 
graphs which remain in the Wess-Zumino model after the chirality condition has been imposed 
and the D-algebra performed yield integrals whose e-expansion involve irrationals or non-zeta 



transcendentals as occurs in the 0(N) (ft theory, [25, 21]. One area where such issues could be 
addressed further is in the calculation of higher order corrections to some of the other critical 
exponents. We would hope to return to this issue later. 
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Fig. 1. Leading order Schwinger Dyson equations. 




Fig. 2. Additional graphs for the a Dyson equation. 
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Fig. 3. Intermediate integral in the calculation of the non-planar graph. 
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Fig. 4. Equivalent three loop integrals. 





Fig. 5. Intermediate integrals in the calculation of the four graph with insertion 
on the central line. 
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